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Abstract 

In conventional beam based alignment (BBA) procedures, the relative 
alignment of a quadrupole to a nearby beam position monitor is determined 
by finding a beam position in the quadrupole at which the closed orbit does 
not change when the quadrupole field is varied. The final focus magnets of the 
interaction regions (IR) of circular colliders often have some specialized prop- 
erties that make it difficult to perform conventional beam based alignment 
procedures. At the HERA interaction points, for example, these properties 
are: (a) The quadrupoles are quite strong and long. Therefore a thin lens 
approximation is quite imprecise, (b) The effects of angular magnet offsets 
become significant, (c) The possibilities to steer the beam are limited as long 
as the alignment is not within specifications, (d) The beam orbit has design 
offsets and design angles with respect to the axis of the low-beta quadrupoles. 
(e) Often quadrupoles do not have a beam position monitor in their vicin- 
ity. Here we present a beam based alignment procedure that determines the 
relative offset of the closed orbit from a quadrupole center without requiring 
large orbit changes or monitors next to the quadrupole. Taking into account 
the alignment angle allows us to reduce the sensitivity to optical errors by one 
to two orders of magnitude. We also show how the BBA measurements of all 
IR quadrupoles can be used to determine the global position of the magnets. 
The sensitivity to errors of this method is evaluated and its applicability to 
HERA is shown. 



1 Introduction 

The new HERA interaction regions are designed to achieve a maximum possible lu- 
minosity by strongly focusing the proton beam. This results in /3-function values at 
the interaction point (IP) which are in the range of the bunch length. This new de- 
sign includes superconducting combined function magnets inside the colliding beam 
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detectors HI and ZEUS which focus the 27.5GeV lepton beam in the vertical plane 
and bend the beam away from the 920GeV proton beam. This allows to place the 
low (3 magnet for the protons as close as 11m to the IP. The synchrotron radiation 
produced by the beam separation has to be absorbed far away from the IP. There- 
fore the vacuum chambers downstream of the IP have a keyhole shape to allow the 
synchrotron radiation fan to propagate through the low-beta quadrupoles. These 
have a 28mm gap between the coils. The aperture of the flat part of the downstream 
vacuum chambers is only 18mm. This is critical because of the height of the syn- 
chrotron radiation that is generated in the upstream low-beta quadrupoles. Due 
to the large vertical divergence of the beam in these quadrupoles, the synchrotron 
radiation fan will only fit inside the keyhole shape if the quadrupoles in the low- 
beta region are aligned to a precision of better than 0.5mm. By optical surveying, 
a precision of about 0.3mm can be achieved under optimum conditions which are 
not given in the interaction region (IR) with shielding walls and a large detector in 
between the two halves of the straight section. Beam-based alignment was proposed 
to verify specifications that cannot be verified to a satisfactory precision by the sur- 
vey procedure. A precision of magnet alignment of 0.1mm appears to be desirable. 
The magnets of the HERA IR are movable via remote control and can be adjusted 
in an iterative way without access to the magnets. 

Beam-based alignment is a technique of deriving the position of a quadrupole 
magnet from the analysis of difference orbits that are generated by the variation 
of the strength of this quadrupole. If the central orbit of the beam is not in the 
center of the quadrupole, the beam experiences a dipole field that changes the orbit. 
Beam-position monitors detect the changes of the orbit around the ring. The offset 
of the beam with respect to the quadrupole axis is then determined by analysis of 
the difference orbit. The result may be used to calibrate the offsets of nearby beam 
position monitors or to mechanically re-align the quadrupole magnets. 

This technique has been invented to optimize the performance of the SLC |], [| 
D, |J. It has also been successfully applied to calibrate the beam position monitors 
in the HERA electron ring, where it was the basis for an orbit steering algorithm 
of minimizing the residual vertical kicks which yielded a record electron spin polar- 
ization ^ |7|]. Future accelerators such as NLC will depend heavily on extensive 
beam-based steering algorithms [§, (J. 

The application of beam-based alignment techniques to adjust the magnet posi- 
tions in the new HERA interaction region however encountered a number of difficul- 
ties and problems. The analysis and the solutions of these problems can be helpful 
for future application of beam-based alignment, especially in interaction regions. 

The difficulties we encountered have to do with circumstances that might be 
considered typical for an interaction region: (a) The quadrupoles are quite strong 
and long. Therefore the thin lens approximation is quite imprecise, (b) The effects of 
angular magnet offsets become significant, (c) The possibilities to steer the beam are 
limited as long as the alignment is not within specifications, (d) The beam orbit has 
design offsets and design angles with respect to the axis of the low-beta quadrupoles. 
(e) Often quadrupoles do not have a beam position monitor in their vicinity. Under 
these circumstances the results are very sensitive to errors and it turned out to be 
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very difficult to achieve the desired precision of the beam-based alignment of 0.1mm. 
Moreover, since the beam cannot be centered in all the magnets simultaneously, a 
global analysis of the magnet positions becomes necessary which uses the results of 
the beam-based alignment measurements in all the IR quadrupoles. 



2 The HERA Interaction Regions 

In the following we describe the HERA Interaction region to the extent relevant 
for synchrotron radiation background and beam-based alignment of the low-beta 
quadrupoles. 

The proton and lepton beams collide head on in the interaction point. The 
two beams are separated by combined function magnets, which start on both sides 
at 2m from the IP. Because of the strong synchrotron radiation power of together 
approximately 30kW generated in these magnets, the layout is not symmetric. On 
the left side from which the lepton beam enters, there is a 3.2m long relatively low 
field superconducting magnet (GO) of 90mm full aperture which deflects the leptons 
by 3mrad and focuses them in the vertical plane as the first lens of a low-beta triplet. 
On the right side of the IP, these functions are provided by a combination of a short 
(1.3m length) but large full aperture (120mm) superconducting combined function 
magnet (GG) and a normal conducting conventional quadrupole (GI) with a length 
of 1.88m. These innermost magnets are complemented by a horizontally focusing 
quadrupole of type GI and a vertically focusing magnet (GJ) of 1.88m length on 
both sides of the IP. The double-doublet structure for focusing the protons starts 
at 11.2m on each side of the IP with a half quadrupole with septum plate. Table 



2A\ shows the main parameters of the HERA IR quadrupoles and the location of the 
positron design orbit relative to the quadrupole axis in the center of each IR magnet. 
For electron/proton collisions the values are slightly different. Due to spin matching 
requirements it has not been possible to use exactly these design parameters for the 
quadrupoles, and in the routinely used optics files the computed path for injection 
and for luminosity operation differ by up to 0.5mm in some quadrupoles. 



3 Analysis of Difference Orbits 

3.1 Closed Orbit Changes due to a Quadrupole Change 

As described before, beam-based alignment is the analysis of difference orbits that 



are excited by a change in the strength of a quadrupole as illustrated in figure |0 

The difference orbit is related to the offset and the angle of the beam orbit with 
respect to the quadrupole axis. Therefore, we will derive the relationship between 
magnet alignment, closed orbit, strength variations and the difference orbit. Follow- 
ing standard text-book procedure, the closed orbit is written in linear approximation 
as 



x 



» = T sfi (/ - Mo)" 1 d Lfi + 4,o • (3.1) 
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Table 2.1: Design parameters for the HERA IR quadrupoles and the offset of the 
positron design orbit in the center of these quadrupoles. 



Name 


/ 

[ m ] 


s 

[ m ] 


focusing 
direction 


k 

[m- 2 ] 


x off 

[mm] 


X 'off 

[mrad] 


QL16L 


1.033 


-54.881 


y 


-0.112026 








QL14L 


1.033 


-42.930 


X 


0.055143 








GJ8L 


1.88 


-9.172 


y 


-0.132197 


-3.288 


-0.002 


GI7L 


1.88 


-6.965 


X 


0.246518 


2.982 


1.322 


GOL 


3.20 


-3.575 


y 


-0.140664 


-5.568 


-0.993 


GGR 


1.3 


2.625 




0. 


-27.913 


-1.217 


GI6R 


1.88 


4.817 


y 


-0.226504 


-10.172 


0.165 


GI7R 


1.88 


7.218 


X 


0.262005 


-0.409 


-0.358 


GJ8R 


1.88 


9.432 


y 


-0.119000 


-9.130 


-0.504 


QL14R 


1.033 


43.934 


X 


0.048787 








QL16R 


1.033 


54.868 


y 


-0.116307 










Figure 3.1: The orbit position relative to the axis of a quadrupole can be deduced 
from the closed orbit change which is created by a change in the quadrupole's field 
strenght. 



3.1 Closed Orbit Changes due to a Quadrupole Change 
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In this expression, x(s) is the vector with closed orbit x and it's derivative x' at 
some longitudinal position s along the design trajectory, which is chosen as reference 
x(s) = 0. T Sj o is the transport matrix from the position s = to position s. With the 
total circumference of the design curve L the revolution matrix at s is M s = T s+ l, s 
and / is the unity matrix. With the focusing strength k, the curvature k of the 
design trajectory in the dipole fields and Ak from dipole field errors and correction 
coils, the vector d S)So describes the closed orbit distortions along the ring according 
to the inhomogeneous equation of motion 

d 2 

— d StSQ + (k 2 + k)d S)SQ = An with d SQtSo = , d' SQSo = . (3.2) 

We assume that the quadrupole to be aligned, the test quadrupole, is the first ele- 
ment in the lattice. For generality the quadrupole can also have a dipole field compo- 
nent. This is important for HERA, since the magnets GG and GO are quadrupole 
magnets with an additional vertical dipole field. For the beam-based alignment 
procedure only the quadrupole strength k is changed, not the dipole field strength. 
Since the quadrupole magnet has a straight axis, the motion through these fields 
is not correctly described with the map of a combined function magnet, where the 
quadrupole field is evaluated around a curved trajectory. In the coordinate system 
which is aligned along the quadrupoles straight axis, the equation of particle motion 
through the magnetic field B y is given by 

x" = = -kx - k . (3.3) 

V 

Charge and momentum are denoted by q and p. The transport map G(x) which 
transports the phase space vector xq from the beginning of the quadrupole to x e at 
its end is given by x e = G(xo) = Gxo + D with 

( cos(lVk) ^= sm(lVk) \ - = _« / 1 - cos(lVk) \ , ^ 
" \ -Vksm(lVk) cos(lVk) J ' k \ Vksm(lVk) J { ' ' 

When the alignment of the magnet with respect to the design trajectory at its 
entrance is described by a shift to Zq and a slope z' Q , and similarly z e describes the 
alignment of the end of the magnet, then the particle transport is described by 

x e = G(x — Zq) + z e + D . (3.5) 

The closed orbit at the end of the quadrupole is given by the periodicity condition 

x e = G(T Lte (x e )) = GiT L . e x e + d L . e - Zq) + z e + D . (3.6) 

After the quadrupole strength has been changed to k + AA; we denote changed 
quantities by a superscript +, and the closed orbit is similarly given by 

f+ = G + (f L , e (x+)) = G + (T Lie xt + 4e - 30 + z e + D + . (3.7) 
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We are looking for a relation between the distance of the closed orbit from the 
quadrupole's center Ax — x — z and the closed orbit change 5x = x + — x. For this 
purpose we use equation ( |3.6| ) to eliminate Ti^Xe + d>L ie — <?o m equation ( |3.7| ) and 
obtain 

x£ = 5x e + Ax e + z e (3.8) 
= G + (T Lie x e + d L , e - z + T Lie 5x e ) + z e + D + 
= G + (G~ 1 (Ax e -D)+T L ,Jx e ) + z e + D + . 

This equation can now be solved to express Sx in terms of Ax, 

(G + G^ 1 - I)Ax e = (/ - G + T Lt£ )5x e - D + + G + G~ l D (3.9) 

This expression can already be used to determine the magnet alignment Ax e at 
its end. However, when neglecting second orders in Ak, the expression become 
simplest when the alignment Ax m in the middle of the magnet is computed. For 
this we introduce the matrix g and vector d which are G and D in equation (|3.-1| ) 
for half the quadrupole length. A useful matrix will be 

(q 1 iVk— sin(iv / fc) \ 

^ +sin( ^) +0(AP), (3.10) 

which reduces the effect of Ak to the center of the quadrupole. In the following we 
will use the abbreviations 

, l^/k + sinly/k _ iVk — sinlVk ...( \o~ \ /011X 

For a defocusing quadrupole (k < 0) sin changes to sinh due to the imaginary unit 
in \J — k, 

_ 1 A(-k) f ^[/vW-rinhavW] , f312 x 

" 2 VW\ \iVW\ + ^HiVW\) o 



With the revolution matrix M m = gTi >e g for the middle of the test magnet equation 
(|3.9| ) leads to 

5cj- l Ax e = (/-(/ + 5)M m )^- 1 5f e - g-\D + - D) + Sg^D . (3.13) 

Relating Ax e to the center of the magnet leads to Ax e = gAx m + d and similarly 
D = gd + d. To leading order in Ak the difference orbit around the ring 8x(s) = 
T s ^ m g~ 1 Sx e is then given by 



Sx(s) = T s , m (J - M m )- l 5ix m -z m - d + AkS^g-'dkD) . (3.14) 



3.1 Closed Orbit Changes due to a Quadrupole Change 
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This shows that for a quadrupole with an additional dipole field, the closed orbit 
distortion Sx is not created by the distance x—z between closed orbit and quadrupole 
axis, but by the distance between the closed orbit and an axis that is shifted by / 
from the quadrupole axis, with 



f = d- AkS-ig-^D = j with / = - I ^-^r 1 -11. (3-15) 

To simplify the notation, we will now use Ax m = x m — z m — / for the closed orbit 
with respect to the modified axis of the quadrupole. The shift of the alignment axis 
by / amounts to — 401/xm for the GO magnet and to — 132/im. for the GG magnet. 

Since T s , m (I — M m ) _1 is the closed orbit generator, the difference orbit Sx(s) is 
created by an effective kick in the center of the test magnet. However, there is not 
only an angle kick 6 m like in the thin lens model of a quadrupole, but there is also a 
position kick A m so that the difference orbit which is created by the change of the 
test quadrupole has two terms, 

Sx(s) = 5xe(s)9 m + Sx A (s)A m , (3.16) 

Sx e = T a , mll {I-M m )^ + T a , ml2 {I-M m )£ , (3.17) 

Sx A = T 8imll (I-M m )^ + T 8<ml2 (I-M m )^ , (3.18) 

6 m = -Akla + Ax m , (3.19) 

A m = Akl^-Ax' m . (3.20) 
k 

Contributions from angular offsets become important if s ™^ffi is significantly smaller 
than unity, which is the case for the HERA low-beta quadrupoles as shown in table 



4.3. The well-known formulas 



T s>m n = W^-[cos(0(s) - m ) + a m sin(0(s) - <p m )} , (3.21) 

y Pm 

T s ,mi2 = \/p(s)P m sm((j)(s) - (f) m ) , (3.22) 
1 / 1 — cos 1i\v + a m sin 1-nv j3 m sin 2tiv 



(J-MJ- 1 
lead to 



4sin 2 7rz/\ — 7 m sin27rz/ 1 — cos 2txv — a m sin 2ttv 



r / \ I qi „\ q cos(|0(s) -0 m | - ?ry) . , 

6X {S) = V/5( S )Pm TT^l \ ' ( 3 - 23 ) 

2 sm \txv) 



x ,s f3(s)a m cos( \(f)(s) -(f) m \-ttu) - sin 0(s) - m \-irv) 
ox A {s) = \ — -— — r .3.24 

V Pm 2sm(7TZ/) 

The contribution Sxe is the conventional closed orbit for a correction coil at the center 
of the test magnet, where (ft = <p m . The contribution Sx A can be compensated by a 
correction coil at d> a = d> m — atan(-M since 

ox A {s) = VP( s )7m — \ sign(a m ) (3.25) 

2 sin txu ) 



s 
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when is not between <p m and <p a . For <p a we take the branch where the atan 
function is in [— |, |]. A closed orbit correction program with correctors at these two 

phases will readily determin 9 m and A m y/^sign(a m ) as proposed corrector kicks. 

These lead immidiately to Ax m and Ax' m with equation ( |3.19 ) and ( 3.20| ). 



3.2 Kick Compensation Method 

Quadrupole errors around the machine might lead to a misinterpretation of the 
quadrupole offsets to be evaluated. We therefore propose to create a closed bump by 
changing the strength of the test quadrupole by appropriately exciting two corrector 



coils as shown in figure p^2 




Figure 3.2: The orbit position relative to the axis of a quadrupole can be de- 
duced from the angles required to close a bump which is exited by changing the 
quadrupole's field strength. 



The difference orbit is thus a closed bump, which starts with 5x — 0, 5x' = 
at the test quadrupole. The amplitude and slope 5x(s) within this bump is derived 
from 

5x(s) = T^G' 1 [(G + - G) (x - z ) + D + — D] , (3.26) 

with the original closed orbit x and the quadrupole offset and angle alignment z 
at the beginning of the test quadrupole. We again refer to the alignment in the 
quadrupole's center by 

5x(s) = T s , m g~ l [g5 3 g~ 1 (x m -z m -d) + D + -D] (3.27) 
= T S)m 5(x m - z m - /) , f = d- Ak5~ 1 g~ 1 d k D . 

The difference orbit vanishes after the second correction coil so that 

^ J = T S2tm SAx m + T S2iSl ( ^ J + ( ^ 2 J ( 3 - 28 ) 

vanishes with x 2 = and x' 2 = 0. Here the deviation of the closed orbit from the 
modified quadrupole center Ax m = x m — z m — f has been used again. The closed 



3.2 Kick Compensation Method 



9 



orbit inside the quadrupole is then 



Ax m = -r 1 ^ ( ° ) +T m , S2 f fl ° )] . (3.29) 



7 2 



To simplify notations, we again use a + = [iVk + sm(lVk)]/(2ly/k), a~ = [iVk - 
sm (lVk)]/(2lVk) and 0im = 0i — 0m- The total transformation between the com- 
pensating kicks = (0i, ©2) and the test quadrupole offset vector is then 



Xrn-Z m = -( ' )+A- 1 e, (3.30) 




A' 1 = —I V ^ 



Where the effective center shift / is given in equation ( |3.15| ) and appears whenever 
the quadrupole field which is changed for beam based alignment is superimposed by 
a dipole field. 

Similarly the corrector angles can be determined from the quadrupole alignment 
by the inverse equation, 




(3.31) 



■-^==(cOS 02m + Ctm SHI 02m} 



k 

^== (cos 0i m + a m sin 4> lm ) \ 



These formulas are accurate up to leading order in The program MAD was 
used to simulate the closed orbit and the kick compensation version of beam-based 
alignment for the HERA IR magnets. The inaccuracy of the reconstructed closed 
orbit deviation due to the neglected higher orders in Ak was shown for all the IR 
magnets to be better than 1.3% for Ak/k < 5%. And it was shown that only second 
order terms in Ak/k contribute noticeably to this small error. 

For the quadrupoles QR16L, QR14L, GOL and GOR the error of the linearization 
is shown in figure |3.3j . The deviation between the alignment x m — z m and the first 
order result (x m — z m )x of formula ( |3.30| ) is plotted against Ak/k on a logarithmic 
scale. The simulations were performed for the displayed range of Ak/k. For even 
smaller Ak/k numerical inaccuracies dominate the compution. The error increases 
linearly, which shows that only next to leading order effects contribute noticeably 
to the errors and that these are small. 

The nature of this transformation becomes more transparent if the two kicks 
9 = {61,62) are replaced by two Active kicks 6f which occur at a betatron phase 
difference of exactly ir and — atan(l/a m ) from the test quadrupole respectively and 
are normalized to the /^-function. We again take the branch where the atan function 
is in [— |, |]. With 0i m — > n and 02 m — > — atan(l/a m ) this leads to 
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Table 3.2: Twiss parameters at positron injection in the center of the quadrupoles 
around the ZEUS IR and in the correction coils right and left of the IR which were 
used to create the closed bumps needed for the beam-based alignment procedure. 
(The strength of the IR quadrupole in the currently used injection optics deviate 
from the design values of table [O] by up to 0.7%.) 

Name j3 x a x <j) x f3 y a y <j) y 



[2tt] [ m ] [2tt] 



CH101L 


14 


.609 





,141 


12. 


,550 


14, 


773 


-1 


,571 


11, 


650 


CV81L 


13 


.568 


2 


,193 


12, 


732 


8, 


,052 


-0 


,260 


11, 


846 


CH75L 


14 


.082 


-0 


,608 


12, 


823 


4, 


426 


-0 


,396 


12, 


,055 


CV56L 


26 


.906 


2 


,734 


12, 


,945 


78, 


231 


-6 


,393 


12, 


223 


QL16L 


23 


.463 


1 


,140 


12. 


950 


86, 


,197 


-1 


,784 


12. 


,225 


QL14L 


33 


.662 





,214 


13. 


,021 


30, 


123 





,841 


12, 


262 


GJ8L 


38 


.355 


-5 


,776 


13. 


,297 


61, 


,265 


6, 


790 


12. 


396 


GI7L 


78 


.538 


1 


,947 


13, 


303 


23, 


,525 


2, 


619 


12, 


406 


GOL 


10 


.899 


3 


,760 


13. 


322 


27, 


998 


3, 


462 


12. 


426 


GGR 


5 


.040 


-1 


,105 


13. 


,614 


19, 


,374 


-7 


,296 


12, 


890 


GI6R 


14 


,276 


-4 


,894 


13, 


,659 


54, 


196 


-1 


,265 


12, 


900 


GI7R 


62 


,721 


-3 


,665 


13, 


,671 


30, 


056 


-0 


,211 


12, 


910 


GJ8R 


33 


,574 


4 


,327 


13, 


,679 


60, 


368 


-4 


,849 


12, 


919 


QL14R 


37 


,975 





,012 


13, 


,897 


17, 


,258 


-1 


,076 


13, 


190 


QL16R 


23 


,143 


-1 


,007 


13, 


964 


78, 


,119 





,934 


13, 


242 


CV56R 


26 


,398 


-2 


,619 


13, 


969 


71. 


,756 


5, 


,345 


13, 


,244 


CH75R 


22 


,662 





,885 


14, 


066 


8, 


,679 


0, 


509 


13, 


369 


CV81R 


19 


,151 


-2 


,690 


14, 


123 


6, 


,176 


0, 


,778 


13, 


,514 


CH101R 


9 


.592 


-0 


,461 


14. 


336 


17, 


,520 


1, 


,879 


13, 


,812 



3.2 Kick Compensation Method 
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for QL16L for QL14L 





Figure 3.3: Deviation between the alignment x m — z m and the result (x m — z m )\ of 
linearization in Ak/k on a logarithmic scale, i.e. log 10 ( ^ m ^ m - >1 — 1) for the accuracy 
of position reconstuction (lower curves) and log 10 (^ 2 f^r^ i — 1) for the accuracy of 
the angular reconstruction (upper curves). 



12 4 INFLUENCE OF MEASUREMENT ERRORS AND IMPERFECTIONS 



— * 

The Active angles Of = (0/i,#/2) are now clearly related to the alignment. The 
kick 6fi corrects the oscillation from the magnet offset so that 0/i = —\f^fi m with 
equation fl3.19| ). The kick and 8f2 corrects the oscillation due to the angle of the 
closed orbit relative to the magnet axis and 9f 2 = v /7mA m due to equation (|3.20|) . 
The relation between the real and the Active angles turns out to be 



/A (c os 4>i, m + Oi m sin 1>m ) a/^(cos0; 
/~{3\\J 1 + a4 sin <f> 1>m 



+ a m sin (p 2 
%yj l + a^ sin (p 2 ,m 




(3.33) 



4 Influence of Measurement Errors and Imperfec- 
tions 



In the following, we will consider the errors in case of the compensating kick method. 



4.1 Error of Thin Lens Treatment 

Let us first consider the error that is made by treating the test quadrupole as a thin 
lens. The effective kick in the quadrupole of equation ( [3.10D becomes in thin lens 
approximation 



leading to — A/c/(x^ m — z m ) = a + (x m — z m — /). Comparing this to equation ( |3.19j ) 



leads to the error of the thin lens version of beam-based alignment measurements, 

x^ -z m = a + (x m -z m -f) . (4.35) 

The error has two components, a scaling error of a~ = 1 — a + which is shown in 
the 2nd column of table |4l| and an absolute error of —fu + which is 452/xm for the 
GOL and 132/mi for the GGR magnet. 

For a HERA IR quadrupole GO the scaling error amounts to 11%. Given the 
systematic horizontal offset of —5.5mm in this magnet and an additional closed orbit 
deviation of up to 5mm, the absolute error due to thin lens analysis could be in the 
order of 1mm. Also for the magnet GI in the HERA IR a thin lens evaluation could 
lead to an error of up to 1mm since the scaling error amounts to 6% and the offset 
could be 10mm design offset plus a closed orbit deviation of 5mm. The error in case 
of a standard lattice quadrupole with k pa 0.1m -2 and I = lm however is rather 
small. Even with an orbit offset of 5mm, the measurement error for such an element 
with zero design offset is only 50/mi. 

The influence of the angle error in the quadrupole alignment is completely ignored 
in the thin lens model. To estimate the relative importance of the oscillation excited 
by the angle alignment, we investigate the Cour ant- Snyder invariant e^ x i of the part 
of the difference orbit which is due to the angle error and e& x which is the part due 
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Table 4.3: Measures of the accuracy of a thin kick approximation for HERA-e's 
interaction region quadrupoles (left) and of the relevance of angular alignment versus 
position alignment given by equation ( f4.36| ) . 

Name o~ = l ^^ sin ( 1 ^) i iVk-smiVk 

horiz. vert. 



QL16L 





010 





006 


mm / mrad 





002 


mm/mrad 


QL14L 





005 





003 


mm/mrad 





004 


mm/mrad 


GJ8L 





038 





044 


mm/mrad 





032 


mm/mrad 


GI7L 





069 





008 


mm / mrad 





036 


mm/mrad 


GOL 





112 





290 


mm / mrad 





104 


mm/mrad 


GGR 










042 


mm/mrad 





054 


mm/mrad 


GI6R 





064 





100 


mm/mrad 





009 


mm/mrad 


GI7R 





073 





018 


mm/mrad 





010 


mm/mrad 


GJ8R 





034 





038 


mm/mrad 





024 


mm/mrad 


QL14R 





004 





002 


mm/mrad 





008 


mm/mrad 


QL16R 





010 





005 


mm/mrad 





002 


mm/mrad 



to the quadrupole shift. With A m and 6 m from equation ( |3.19| ) and ( |3.20| ) we obtain 



^Ax' Im^m = 1 Q_ V 1 + °4 Ax' 

6Ax V Pm9 2 m k CT+ (3 m Ax 1 ' J 



For the HERA IR magnets the 3rd column of table |4.3| shows this ratio of oscillation 
amplitudes. 

For the long, superconducting GO magnets in the HERA IRs the oscillation 
amplitude e& x / is 29% of the oscillation amplitude due to the quadrupole shift. The 
quadrupole's angle and the corresponding error from neglecting it are therefore quite 
significant. Thus even if the contribution of angular alignment is considered small, 
one should take into account only that component of the difference orbit which has 
the proper phase relation to the test magnet in order to avoid large errors. This 
is especially important, if there is a large value of a m in the center of the test 
quadrupole. 

4.2 Influence of Optical Errors 

Beam optics distortions between the compensating kicks in a beam-based align- 
ment measurement and the test quadrupole lead to misinterpretation of the differ- 
ence orbit and a corresponding error of the evaluation. Given optical errors 8(3 = 
(5p m , 5 j3i, 8(52), 8a m and S(j) = (4>i m , 4>2m), the change Ak in the test quadrupole re- 
quires correction kicks 9 to close the bumps which are obtained by inserting the per- 
turbed optical functions in to equation ( |3.31| ). The inferred orbit in the quadrupole 
is however obtained by equation ( [3.30D with the unperturbed optical functions and is 



therefore erroneous. This erroneous result of the beam-based alignment procedure 
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Table 4.4: Errors of the beam-based alignment procedure for the HERA IR 
quadrupoles at injection due to a phase error of 5<p = 0.01 ■ 2tt within the closed 
bump of the kick compensation method. 



in ame 


O A x LAX — 1 


Pt . , A T err 
OAx'^X 




Pi. .At'™ 1 






mm / mi au 


mi ao / mm 






-U.U ( 


n nn 
u.uu 




n ri7 

U.U ( 


OL14L 


-0.01 


0.00 


-23.74 


0.01 


GJ8L 


0.36 


0.02 


-8.312 


-0.36 


GI7L 


-0.12 


0.00 


-16.29 


0.12 


GOL 


-0.24 


0.07 


-0.843 


0.24 


GGR 


0.07 


0.00 


-2.249 


-0.07 


GI6R 


0.31 


0.03 


-3.129 


-0.31 


GI7R 


0.23 


0.00 


-12.99 


-0.23 


GJ8R 


-0.27 


0.01 


-7.264 


0.27 


QL14R 


-0.00 


0.00 


-26.79 


0.00 


QL16R 


0.06 


0.00 


-16.42 


-0.06 



is here referred to as Ax 011 . We refer to the matrix in equation ( |3.31| ) as A(@, a r 
The matrix in equation ( |3.30|) is A -1 , 



AxZ T = A-\(3, a m , $)A0 + 6(3, a m + 6a m J+ 5<f)Ax . (4.37) 

For simplicity let us now assume a phase error, so that the a and /3 functions 
do not change and we assume 6<j) = £0i,m = 64>2,m which means that no optics error 
occurs between the corrector magnets. Here we will neglect all nonlinear terms in S(p 
by replacing cos(<50) by 1 and sin(<50) by 6(f>. With the equations ( |3.31| ) and ( |3.3(J| ) 
the result can be expressed in the following way 

AxZ r - Ax m = 6<p( ~f m CT+ l * lm °- ) Ax m . (4.38) 

The term in the position error which is proportional to Ax m is thus simply given by 

d Ax Ax™ - 1 = -OlrM . (4.39) 

With a phase deviation A<p = 0.01 ■ 2n this error is 24% for the GO quadrupole. For 



the same phase deviation, table [4.4| shows all these errors for the HERA IR magnets. 
The error in the position measurement that is introduced by the angle alignment is a 
few percent. However, the term that generates the error in the angle determination 



is shown to be huge in the third column of table |4.4j This will prevent a precise 
measurement of the angular alignment. 

The here studied case of a pure phase error is somewhat artificial. When the other 
optical functions are also perturbed, then the evaluation becomes rather elaborate 
and the errors depend strongly on the location of the optical element that courses 
them. We now assume that there is one thin lens quadrupole error with focal 
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strength Ski at position q in between the test magnet and the two correction coils. 
The kicks 9 in the correction coils are then related to the alignment by equation 

Ax m = -5T l T m , q f _} 5h \ \ [T q , sl T q , S2 [I)]- ( 4 - 4 °) 

The erroneously determined alignment Aa; err does not take the optical error into 
account, 

Ax% = -S^TmAT^ ( £ ) + T ^ ( 2 )l ■ ( 4 - 41 ) 
We therefore obtain the relation 

Ax^ - Ax m = 5T l T mA ( ^ ° J T ? , m 6A£ m . (4.42) 

The error is a linear combination of the deviation Ax m from the magnet center and 
the deviation of the slope, 

(Ax™ - Ax m = Ax m {d Ax Ax™ - 1) + Ax' m d Ax ,AxZ r • (4.43) 

and similarly for the error of the angular alignment determination. The exact value 
of the terms in the matrix that relates Ax° n and Ax depend on the optical param- 
eters, especially on the phase advance between the error and the test quadrupole. 
When one inserts as a worst case scenario for each of the matrix elements the phase 
4™ where it has the maximum absolute value, one obtains 



1 



Max\d Ax Ax e ™- 1| = IShlPt-daJ + y/l + a 
Max|^Ax e ^| = l^l/Vy. 



in > 



m ka~ 



krr + 



a' 



Max|9 A:[ ,Ax^ r -l| = \5h\f3 q ^(\a m \ + y/l + o&) . (AAA) 



These values are shown for the HERA IR in table |4.4j when a thin lens quadrupole 
error with tune change ^bk\f5 q = 0.01 is assumed. Table fO| shows which relative 
quadrupole errors in the IR lead to such a tune shift. 

For specific quadrupole errors in the HERA interaction region, the sensitivity of 
this beam based alignment procedure was also evaluated. Table [1.7| shows the error 
of the constructed beam offset in the GJ8L, GI7L, and GOL magnet which occurs 
when the field strength in one of the other quadrupoles on the left side of the IP has 
an error which leads to a tune shift of 0.01. We only show the most relevant term 
of the error, d Ax (Ax cn — Ax). 
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Table 4.5: Maximum of the errors of the beam-based alignment procedure for the 
HERA IR quadrupoles at injection due to some focusing error at position s q with 
tune shift -^f3 q 5ki = 0.01. Corresponding quadrupole errors are shown in table [O. 



in ame 


\Ft* A T crr 1 I 
\Oa x /ax — l\ 


\Of\ x 'L\X | 


\rt. Ai»' crr l 
\0/\ x lAX 








111111/ 1111 clU. 


1111 cLU. / 111111 




OL16L 


0.17 


0.00 


33 


0.17 


QL14L 


0.08 


0.00 


47 


0.08 


GJ8L 


0.73 


0.03 


17 


0.73 


GI7L 


0.26 


0.00 


33 


0.26 


GOL 


0.48 


0.14 


2. 


0.48 


GGR 


0.16 


0.01 


4. 


0.16 


GI6R 


0.62 


0.06 


6. 


0.62 


GI7R 


0.47 


0.01 


25 


0.47 


GJ8R 


0.55 


0.02 


14 


0.55 


QL14R 


0.06 


0.00 


53 


0.06 


QL16R 


0.15 


0.00 


32 


0.15 



Table 4.6: Relative filed strength errors for the IR quadrupoles which lead to a tune 

shift AQ X or AQ y of 0.01. 

Name QL16L QL14L GJ8L GI7L GOL~ 

4^6% 6^5% L3% 0Wo 2.6% 

1.3% 7.3% 0.8% 1.1% 1.0% 

Name GI6R GI7R GJ8R QL14R QL16R~ 

2l% 04% L7% 6^% 4.5% 

0.5% 0.8% 0.9% 14% 1.3% 



Table 4.7: Horizontal errors d& x x eTT — 1 of the beam-based alignment procedure for 
three HERA IR quadrupoles at injection due to a error of the field strength in one 
of the other IR quadrupoles which leads to a 0.01 tune shift. All magnets were 
assumed to be correctly aligned. 



error element 


for GJ8L 


for GI7L 


for GOL 


QL16L 


-0.432 


0.217 


0.309 


QL14L 


-0.687 


0.259 


0.463 


GJ8L 




-0.005 


-0.008 


GI7L 






-0.008 
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4.3 Reduction of Sensitivity to Errors 

One source of errors is an imperfect determination of the compensation kicks 6\ and 
02- Since the determination of Ax' m is very prone to errors, as can be seen in the 
tables fO and [4.5|, it is not worth trying to determine the angle alignment. But 



we will assume that the angle alignment Ax' m = x' m — z' m is approximately correct 
and we will therefore require our compensation to lead to the design value of the 
angular alignment. While the angles are measured, we assume that the correct 
angles to close the bump would have been — A0. Since the errors A0 are not 

— * 

known, we introduce an estimate A0* of the erroneous angle such that equation 
( p.30| ) leads to an estimated alignment of 

A** m = A~ l 0, a m , 0)[0 - AO*] = ( A fj \ . (4.45) 

With a 2 = (A^^A^l) we write a 2 • [0 — A©*] = — z'£. This condition should be 

satisfied for a set A0* of angles which is as small as possible, i.e. |A0*| 2 should be 
minimal. We can use Lagrange multipliers to minimize, 

AO* 2 + A[a 2 • (0 - AO*) + z°] -> minimum , (4.46) 
2A0* - Xa 2 = , (4.47) 
a 2 • AO* = a 2 • 6 + z'° . (4.48) 

These equations lead to AO* = a2(a2-0 + ^)/|a2| 2 . With di = (A^l, A~{ 2 ) equation 
Q3.30Q determins the alignment to Ax m — a± ■ 0. When the above estimate AO* is 
used, the estimated alignment is given by 

1 ri — * 1 2 —*V t — * j v _ fp-\ s~t CL\ ' Ct^ 1 q 

Ax m = -—[a 2 a 1 -{a 1 -a 2 )a 2 Q--—-z m 
\a>2\ l a 2T 

-a 2 (a 2 , -a x ) J, - -r=-r^z m 



\a 2 f \ a 2 J |cta 

1 7?T ( A 2~l ~ A ll \ ( A ll A 12 \ 3 «1-«2 J0 
^"2 /l-l A-l A-l A-l U 1-19 Z n 



\7J 1 2 "2 I A-l _ A-l j 1 A-l A-l I I- 12 ^rn 

\a>2\ \ ^22 A 12 J \ ^21 ^22 / \0-2\ 

det(A~ 1 ) ^ T f -1 \ * di-a 2 > 



(4.50) 



With equation ( |3.30| ) for A 1 this yields 



Ax * _ VA&sinfoi y^sin^mQi - y/Jhsiri(j) lm 2 
Akla + a/A^ Pi sin 2 <f> lm + p 2 sin 2 <p 2m 
>o cr- p x sin (f) lm cos im + p 2 sin (p 2m cos 2m 

kfj m a + fix sin lm + sin 2m 

When the determination of the angles has an error with standard deviation er@, 
then the errors in the determination of Ax m and Ax' have the standard deviations 
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Table 4.8: Systematic error Ax*^ — Ax^ of the procedure which reduces the mea- 
surements sensitivity to the kick angles 0. This error increases linearly with the 
deviation from the orbits design angle in the magnet. 



GJ8L 


GI7L 


GOL 


GGR 


GI6R 


GI7R 


GJ8R 


m.rad 


q ion 
m.rad 


291 

m.rad 


rnrad 


-99 ^ m , 

m.rad 


rnrad 


37 ,m , 

rnrad 



|ai|o"e and (a^cxe when equation ( |3.30|) is used, leading to 




COS ( 



'lm 



a m sin0i m ) 2 + 



fa 

Prr, 



cos</> 2m + a m sin, 



'2ml 



(4.52) 



When the equation ( |4.49j ) is used, the standard deviation of Ax* m is always smaller, 

det(A-^/0 -1\ 



0-q 



\a 2 \ 



1 ->T 

-a 2 



-1 

1 



a//5i/3 2 sin 02i 



Akla+^, ^ S i n 2 lm + ^ 

sin^ 



^2m 



Especially for large a m his spread of results in Ax* m is drastically smaller than the 
spread in equation ( [4.521) . If the angular alignment of the orbit relative to the 



magnet is not the design value 



Jo 



then Ax* contains a systematic error since 



A© = does not lead to the correct alignment Ax r 
it leads to 

A^ m — Ax m — (Ax m + z m ) ; — ; - 



Ax n , 
a' 



ka + j3„ 



[air 



\a>2\ 

Z m "T" z m) X 



Si ■ but with a 2 • © = Ax' 



(4.54) 



/?i cos 0i m sin 0i m - f3 2 cos 2m sin <f) 2m - 



/?i sin 2 lm + (3 2 sin 2 



(4.55) 



^2m 



This systematic deviation of Ax*^ = Ax^(A0 = 0) from Ax Q m = Ax m (AQ = 0) is 
shown in tabel L8 for an angular deviation x' of 1 rnrad. Figure | | shows that a 



spherical error distribution for Gi and 02 leads to an elliptical distribution for Ax m 
and Ax' m . The large spread in Ax m is reduced by the estimation of A0* as shown 
in the figure. This, however, introduces the systematic error x*^ — x® m which is also 
shown. 

Now we will show that this procedure also reduces the sensitivity to optical 
errors. In the tables |4.4j , |4.5| and |4.7| is has been seen that the most important error 
of the alignment determination is due to the term d/^ x Ax c " — 1. This is mostly 
due to the fact that the Twiss parameter a m that contributes to this term can be 
relatively large. We will now show that the here proposed method of error reduction 
makes this term independent of a for all types of optical errors. For an alignment 

Ax m , the angles that close the bump are given by = A(j3, a m , <p)Ax m where 
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the tilde indicates Twiss parameters which are perturbed due to an optical error. 
With equation ( 4.53Q the estimate of the alignment is computed by 



detjA" 1 ] 



ai • a2 'o 

7, 

I -* 19 m 

\a 2 r 



The most disturbing error contribution d Ax Ax^ — 1 is then given by 



(4.56) 



d Ax AxZ r - 1 = 



det(A- 1 ) 



fl 2 



— — z — Sin <p2m 

■—— — sin <pi m 




(4.57) 



Since det(A _1 ) as well as | c? 2 1 2 do not depend on a m , the error contribution no longer 
depends on a m , no matter which optical perturbation occurs. Equation ( |4.38| ) and 
table [14| show the error terms introduced by an optical error that only changes the 
betatron phase. When the error reduction method is used, the error of the alignment 
determination can be computed from equation ( |4.57| ) to be 



d Ax Ax 



err 

m 



(3\ sin 2<pi m + j3 2 si n 20 



2m 



f3\ sin^ 



+ p 2 sin^ 



(4.58) 



■>2m 



Also the term d Ax Ax 
For <p lm = ix/2 and 2m 



err 
m 



1 depends on the Twiss parameters at the corrector coils. 
: 7r we obtain 



Aa;* 



-Axi 



(4.59) 



The error becomes completely independent of Ax m and the already small error due 
to Ax' m in equation ( |4.38| ) is reduced by \A + a m - For the phase advances which 
are realized in the HERA IR between the test magnet and two horizontal correction 
coils at 101m and 75m left of the IP, the error terms are shown in table [4.9| . The 
error has been reduced to less than 1% for all magnets, whereas it was up to 36% 
without error reduction. 

The maximum error terms that can occur due to a focusing error somewhere in 
the bump depend also on the phase advances when the error reduction method is 
used. For <p lm = n/2 and 2m = n one obtains 



Mtxx\d Ax Ax* m - 1| 
Max|d A:r ,A<J 



6^(1 +yr+ ai) g- 

2 k(3 m a- 



(4.60) 
(4.61) 



Both terms are always smaller than the maximum errors in equation (]4 Ej) without 



error reduction. Especially the first error term is significantly smaller as can be 
seen in table |4.9| where these maximum errors are plotted for the IR of HERA. In 
table [4. 1Q| the error for one of the IR magnets is shown which occurs when another 
IR magnet causes the focusing error. With error reduction also these errors are 
significantly smaller than those in table |T7. 
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Table 4.9: After error reduction: the horizontal beam-based alignment procedure 
for the HERA IR quadrupoles at injection due to a phase error of 5<p = 0.01 ■ 2ir 
within the closed bump of the kick compensation method. 



Name 


d Ax Ax err - 1 


a Ax /Ax crr 

mm/mrad 


Max\d Ax Ax err - 1| 


Max\d Ax >Ax ei1 
mm / mrad 


QL16L 


-0.002 


0.000 


0.063 


0.000 


QL14L 


-0.005 


0.000 


0.063 


0.000 


GJ8L 


0.009 


0.001 


0.063 


0.003 


GI7L 


0.009 


0.000 


0.063 


0.001 


GOL 


0.009 


0.002 


0.063 


0.023 


GGR 


-0.007 


0.002 


0.063 


0.004 


GI6R 


-0.003 


0.001 


0.063 


0.008 


GI7R 


-0.002 


0.003 


0.063 


0.001 


GJ8R 


-0.001 


0.001 


0.063 


0.003 


QL14R 


0.005 


0.000 


0.063 


0.000 


QL16R 


-0.004 


0.000 


0.063 


0.000 



Table 4.10: After error reduction: d^xAx " — 1 for three HERA IR quadrupoles 
at injection due to a error of the field strength in one of the other IR quadrupoles 
which leads to a 0.01 tune shift. 



error element 


for GJ8L 


for GI7L 


for GOL 


QL16L 


0.048 


0.049 


0.053 


QL14L 


0.004 


0.008 


0.021 


GJ8L 




-0.005 


-0.020 


GI7L 






-0.015 
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We conclude the error considerations by realizing that there are no chances to 
determine the angular alignment of the quadrupoles with the desired precision of 
Ax' m ~ 100/urad. The large sensitivity to optical errors and to corrector settings of 
the estimated quadrupole offset can be reduced drastically by one to two orders of 
magnitude. 

5 The Global Positions of Magnets 

5.1 Combining BBA Data of all IR Magnets 

Since it is not possible to steer the beams to the middle of all quadrupoles for a 
misaligned interaction region, the measurement of the position and angle of the 
beam with respect to a single quadrupole magnet does not give enough information 
to determine the global alignment of this magnet. All the quadrupole offsets and 
angles with respect to the beam have to be determined, and the beam orbit has to 
be consistently modeled, thereby fixing the absolute magnet positions. In order to 
achieve that, the following procedure has been established: the beam offsets with 
respect to all the quadrupoles in the IR are measured for two or more different 
quadrupole settings in the IR. Then a model of the IR that has the initial orbit 
values at the entrance of the IR and the magnet positions as free parameters is 
fitted to the set of measurements. An additional constraint in the fit is that the 
magnet position deviations from their nominal values should be minimal in order to 
connect to the machine coordinate system avoiding a global, unrealistic offset. 

In order to perform this task, we need an explicit formula for the beam orbit as 
a function of alignments and initial conditions. To arrive at such a formula we write 
the transformation of the beam orbit from the center of a quadrupole to the center 
of its neighbor quadrupole as 

[9n-l{%n-l — ^n-l) + G?n-1 + ^n-l^n-l] + <^n,n-l — ^ d"n ~\~ Z n . 

The 4x4 matrix D n n _i is the transport matrix from the end of the n-lst test 
magnet to the entrance of the nth test magnet. On this distance the closed orbit 
distortions d n) n-\ due to corrector coils or field errors are being accumulated. The 
matrix t n describes a drift with half the length of the nth magnet. It is used to obtain 
the alignment at the end of a magnet as t n -i"?n-i or at the beginning of a magnet 
as t~ x z n . In the following it will use the 4x4 matrix T nn _! = g n D nn _x9n-\- The 
matrix T nj „_i transforms from center to center between two neighbored quadrupoles. 
The vector x n describes the orbit and z n is the vector of magnet alignments, both 
taken in the center of the magnet, g n transforms through half the quadrupole with 
index n. This equation can be simplified by using 5x5 matrixes, where the fifth 
column is used to describe the closed orbit distortions. The orbit vector then has five 
components, (x, x', y, y', 1), the alignment vectors z n have in their 5th component. 
The closed orbit deviations d are then all absorbed in the 5th columns so that, after 
combining the terms, the recursive orbit formula reads 

%n — T n:n _\x n _\ — T nn „ 1 (J — (7 n 2]£ n _i)z n _i + (1 — g n t n )z n . (5.62) 



5.1 Combining BBA Data of all IR Magnets 



23 



This recursive formula leads to the explicit expression 

x n = T nfi x - T n>0 (I - goHojzo (5.63) 



n-l 



Using 



Pn,o = -T nj0 {I - g~H ) , (5.64) 
Pnj = T^igjhj - gjtj 1 ) for < j < n , (5.65) 

Pn,n = -9J- 1 ■ (5.66) 

We finally obtain for the orbit in each of the N test magnets the desired form 

u 

X n — Z n = T n flXo + Pn,jZj ■ (5.67) 

j=0 

On the left side appears the expression which is obtained as the result of the mea- 
surement, the right hand side contains the parameters to be fitted, the magnet 
offsets and angles and the initial orbit coordinates. Since there are more parameters 
than measured values, this expression can only be solved by fitting at least two 
different measurements with different quadrupole settings (index m, thus different 
matrices and P^) simultaneously, or adding as additional constraint that the 
magnet positions should differ as little as possible from their nominal value z n = z\. 
Additional constraints can be the readings of the beam position monitors BPMs in 
the IR region. 

From the previous paragraph, it is clear that the angular alignment cannot be 
determined with satisfactory precision. Since the angles of the magnets with respect 
to the beam need to be taken into account for the fit, the design angles of the magnet 
are used. Therefore we can only make use of the position part of the vector equation 
( {5.67] ). This introduces certain errors into the alignment reconstruction that are 
analyzed in the next section. 

We now define a new relationship between the measured values of x n — z n , the 
magnet offsets z Q , z\ . . . Zn and the initial orbit values xo, x' Q . For this we define 
new vectors. The first one includes several sets (m £ {1 . . . M}) of relative position 
measurements (x m — ^ m )^ m - ) in all the IR magnets, as well as the design position z® 

v = ((x - ^o) (1) , (xi - ^i) (1) , ...(x N - z N )W, ...(x N - z N f M \ zl ...z° N ). (5.68) 

The second vector contains the parameters to be determined 

g = (z , ...z N , x { o\ x (1) , . . . x { M \ x (A/) ) , (5.69) 

and the third vector w contains the parameters which are kept fixed, the design 
angles of the magnets and the effects of the dipole corrector settings of the mth 
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measurement, 



w 



(m) 



pirn) 
-n,0 



3=0 



+ E \ p of ■z' 3 °iovje{l,...N.M}, 

1,5 ^— ' L ,J J 1,2 J 



(m) 



else . 



(5.70) 

Note the indices outside the square brackets denote the matrix element, the indices 
inside the square bracket denote the matrix. The measurements are then related to 
the parameters by 

v = Ag + w. (5.71) 
The matrix A contains the matrix elements which are determined by equation ( |5.67| ), 
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(5.73) 
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The solution of the fit with a quadratic norm is 

g= (A T A)~ 1 A T (v - w) . 

In cases when the inverse cannot be computed due to a bad condition of the matrix 
a singular value decomposition (SVD) can be helpful. 



5.2 Error in Magnet Offset Determination by Ignoring the 
Magnet Angle Offsets 

The quadrupole magnet's angular alignment produces a contribution to the closed 
orbit that is to be reconstructed. Since we are not able to measure the angular 
alignment, it is not contained in our model. Therefore, the missing angular offset 
in the model is compensated by an additional false offset in neighboring magnets. 
The magnitude of this error is estimated in the following. 

The effect of the magnet angle on the beam trajectory is described in equation 
( |5.63| ) for a one of the orbit planes by two dimensional sub-matrices, 
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We use the same notation as before, T n j is the transport matrix between the middle 
of quadruple index j and n, gj is the transformation through half of this quadrupole 
and tj is a transformation though the drift of the same half length, z'a is the angle of 
the beam with respect to the design curve. In our model, the missing effect of the 
angle alignment in equation ( [5.80 ) is produced by additional, false magnet offsets 



of neighboring quadrupoles. Quadrupole i requires an additional shift Azi and an 
additional angle Az[ to produce the effect of the angle of quadrupole j, 



X n = T ni Qi ( ^ (5.81) 



These false offsets are given by 
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3 

For two quadrupoles with a phase distance of <frij = — arctan(l/aj), taking the atan 
function in {0,7r], one obtains 

=T l3 = ^7 



This leads to Az[ = so that the quadrupole i alone can compensate the missing 
angle of quadrupole j. The error of the position reconstruction for this quadrupole 
is given by 



ij l iV^ cos dV^) - sin (| y/fy), i 



V Pi ^kikjSm(^Vki) 

The factor between Azi and z'j — z'® depends very much on the values of a in the 
centers of the two quadrupoles an on the chosen phase advance which allows that a 
single quadrupole can compensate the ignored angular alignment of magnet j. 

In a realistic setting, at least two quadrupole positions, i and q, will be recon- 
structed erroneously to take account of the missing angle alignment of quadrupole 
j. For this the following equation has to be satisfied, 

TniQi ^ q 1 J + T nq Q q ^ o 9 ) = ^ nj Q j (y z '- - z'° J (5.85) 
Solving for the position errors leads to 

Az q ~ {j j) sin, 



(5.86) 



Table |5TTT| shows what error an angular alignment error of HERA's IR quadrupoles 
can have on the estimated position of the two neighboring IR magnets. Even for an 
angular alignment error of lmrad, the errors in the reconstructed position is always 
below 150/im. 
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5 THE GLOBAL POSITIONS OF MAGNETS 



Table 5.11: Effect of an angular alignment error of HERA IR magnets on the esti- 
mated alignment of the two neighboring quadrupoles. Since the nominal quadrupole 
strength of GGR is zero, this magnet is not considered here. 

Name Quad to the left Quad to the right 
mm/mrad mm/mrad 



QL14L 


0.000 


-0.001 


GJ8L 


0.029 


-0.017 


GI7L 


0.074 


-0.056 


GOL 


0.104 


0.050 


GI6R 


-0.016 


-0.041 


GI7R 


0.070 


-0.138 


GJ8R 


0.014 


-0.030 


QL14R 


0.001 


0.000 



5.3 Error Propagation in the Fitting Procedure 

The magnet positions as the result of a fit of the measurements can be written as 

N 

z n = ^Bnji/j . (5.87) 
3=1 

The z n are the first N components of the vector g, the yj are the components of 
v — w and B = (A T A)~ 1 A T is the matrix of the least square fit in equation ( 5.77|) . 
We now introduce sets of random errors of the input variables labeled by a 

N 

z n + Az a n = J2(B nj + ABZj)( Vj + Ay") . (5.88) 

3=1 

which gives the error of the magnet positions (neglecting second order terms) 

N 

Az: = E( Ai ^ • Vi + B ^ ■ A vf) ■ (5-89) 

3=1 

We now calculate the expectation value of Az n by squaring the expression and by 
subsequently averaging over the error set a. We assume that for random errors 
the following correlations hold: < AB^AB^ > a = AB? ms 5 jjk and < AyfAy% > a = 
Vrmsfi j,k- With this we finally obtain 



Az n , rms = J2 ^ ■ AjB -s + E B h ■ A ^ms • (5-90) 



3 3 



The first sum has typically values of 1 to 3 when evaluated for the HERA IR. Thus 
an error of a single measurement A(xj — Zj) = 0.1mm propagates, yielding approx- 
imately an error of 0.3mm in the reconstructed magnet position. However, these 
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are only first observations for the case of HERA where the here presented version of 
beam-based alignment is currently being heavily used in the commissioning process. 
This particular application and experiences with this method will thus be reported 
in a separate paper after the successful commissioning of the HERA luminosity 
upgrade. 

6 Conclusion 

We have introduced a beam based alignment method for a general class of combined 
function magnets that can be encountered in collider interaction regions. While 
this method can in principal determine alignment angles, we have shown that these 
angles would be very prone to measurement errors. We have therefore introduced 
a procedure to use the angular alignment to strongly improve the accuracy of the 
position determination by one to two orders of magnitude. Furthermore a proce- 
dure has been presented to determine global magnet positions after the closed orbit 
deviation from the quadrupole axis has been measured throughout the interaction 
region. 
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